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Anomalous minimum and scaiing behavior of 
iocaiization iength near an isoiated flat band 

L. Ge^’^* 


Using one-dimensional tight-binding lattices and an 
analytical expression based on the Green’s matrix, 
we show that anomalous minimum of the localization 
length near an isolated flat band, found previously for 
evanescent waves in a defect-free photonic crystal 
waveguide, is a generic feature and exists in the An¬ 
derson regime as well, i.e., in the presence of disorder. 
Our finding reveals a scaling behavior of the localiza¬ 
tion length in terms of the disorder strength, as well as 
a summation rule of the inverse localization length in 
terms of the density of states in different bands. Most 
interesting, the latter indicates the possibility of having 
two localization minima inside a band gap, if this band 
gap is formed by two flat bands such as in a double¬ 
sided Lieb lattice. 


1 Introduction 

Systems that exhibit flat bands have attracted consider¬ 
able interest in the past decades, including optical [1,2] 
and photonic lattices [3-7], graphene [8,9], superconduc¬ 
tors [10-13], fractional quantum Hall systems [14-16] and 
exciton-polariton condensates [17,18]. As the name sug¬ 
gests, a flat band is dispersionless, or in other words, its 
density of states (DOS) diverges at a particular energy, 
known as the flat band energy. A quantity strongly influ¬ 
enced by DOS is the the localization length in a disor¬ 
dered system, and it displays quite different properties 
depending on whether a flat band is present. For exam¬ 
ple, previous studies in flat band systems have shown 
inverse Anderson transition [19], localization with uncon¬ 
ventional critical exponents and multi-fractal behavior 
[20], mobility edges with algebraic singularities [21], and 
unusual scaling behaviors [22,23]. 

Despite all these interesting findings, there is one im¬ 
portant aspect that has not been investigated systemati¬ 
cally, i.e., the energy dependence of the localization length 
near an isolated flat band, which is separated from its 


neighboring band on each site by a band gap. To be clear, 
here by localization we refer to the finite spatial extension 
of the wave function. There are two mechanisms typical 
in a periodic system that can lead to wave function con¬ 
finement, one by evanescent waves inside a band gap 
[24], where the DOS (of the propagating states) vanishes. 
The other is Anderson localization [25], where disorder, 
ubiquitous in naturally formed and fabricated materials, 
suppresses wave diffusion. 

One usually finds a minimum of the localization length 
near the center of a band gap (“midgap"), where it is most 
unlikely to couple to the neighboring (propagating) bands 
[26-30]. An exception to this conventional understanding 
was reported in a photonic crystal waveguide [31], formed 
by removing one row of unit cells from a two-dimensional 
(2D) periodic dielectric slab, ft was found that the shortest 
evanescent tail in the transverse direction, perpendicu¬ 
lar to the removed row, is at a frequency shifted signifi¬ 
cantly away from the midgap. This observation was later 
explained using the complex wave vector inside the band 
gap of the perfect 2D slab (before the row removal) [32], 
where a flat band was found on one side of the band gap 
as the result of a Van Hove saddle point singularity. 

In this work we show that these anomalous localiza¬ 
tion minima are not limited to evanescent waves. Instead, 
they are a generic feature of an isolated flat band and inde¬ 
pendence of the origin of the latter. Hence these anoma¬ 
lous minima exist in the Anderson regime as well, i.e., in 
the presence of disorder. We exemplify this finding using 
a one-dimensional (ID) tight-binding Lieb lattice, and we 
attribute these anomalous minima to the compact states 
in the flat band [2,18], which is supported by an analytical 
expression of the inverse localization length. This expres¬ 
sion reveals a scaling behavior of the localization length 
in terms of the disorder strength, as well as a summation 
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Figure 1 Band structure of a simple periodic lattice (a) and 
a Lieb lattice (b). In (a) Eb - G, Eq- -G, and G = 1. In (b) 
Ea = Eb = Ec = 0, G = J = 1, and the flat band is located at 
EpB = 0. a is the lattice constant. 

rule of the inverse localization length using an integral 
transform of the DOS in the flat hand and the opposing 
band. We find the latter term to be almost a constant near 
the flat band energy (denoted by EpB below) for weak dis¬ 
order, and the flat band term shows a non-monotonic 
energy dependence that leads to the anomalous minima 
of the localization length. Most interestingly, we show that 
the localization length can exhibit two minima instead of 
one inside a band gap, if this band gap is formed by two 
flat bands such as in a double-sided Lieb lattice. 


2 Properties of the localization length 

2.1 Evanescent waves in the absence of disorder 

Before discussing the anomalous localization minima in 
the presence of disorder, we first briefly review the proper¬ 
ties of evanescent waves inside a band gap in the absence 
of disorder. A band gap, by definition, does not support 
any extended or propagating states, fn other words, the 
wave vector inside a band gap takes a complex value, and 
its imaginary part gives the inverse localization length of 
the evanescent wave [24]. From this point of view, two 
systems with identical band structures (except at some 
singular values of the energy), should also have identical 
localization properties for evanescent waves inside the 
band gap. 

As a case in point, we compare a simple ID periodic 
lattice and a Lieb lattice in Fig. 1. They have two identical 
dispersive bands given by 

£= ±\/G2-t 2/2(1-r cos fca), (1) 

with the Lieb lattice featuring an additional flat band at 
EpB = 0. These band structures are calculated using a 


tight-binding model, with nearest neighbor couplings / 
and G shown in the insets of Fig. 1. Take the Lieb lattice 
for example, the tight-binding model is given by 

i A.j = Eaj Ay -t- GBj , (2) 

iBj = EbjBj + GAj + J[Gj -1- Gy_i), (3) 

iCj = EpjGj + J(Bj + 5j+i), (4) 

where the overhead dot denotes the time derivative and 
we have taken h = l. Xj is the value of the wave function 
on the jth site of type XiX- A,B,C\j = 1,2,...), and/j^y 
is the energy on site Xj . Clearly the wave function and 
hence the localization length do not depend on the overall 
energy scale, because all the terms in Eqs. (2)-(4) are linear 
in energy. We can thus use an energy scale at will, which 
we choose to be G and set it to be 1. 

The model above assumes noninteracting particles, 
with which the localization of fermions, bosons, and clas¬ 
sical waves can be treated on the same footing in their 
respective regime of validity, i.e., as a general wave phe¬ 
nomenon where propagation is suppressed [33]. We also 
note that although Eqs. (2)-(4) are written in the form of 
the Schrddinger equation, it applies to the wave equation 
as well in the paraxial regime, which can be realized, for 
example, by coupled parallel waveguides [3,4]. There the 
time derive is replaced by the spatial derivative along the 
axial (propagation) direction, and the energy is replaced 
by the propagation constant [34]. 

To find the localization length of the evanescent wave 
inside the band gap(s) shown in Fig. f, we simply need to 
invert Eq. (1) and solve fc as a function of E: 

k(E) = 7i±ln\a[E) + V(r^iE)-iy a{E] = l+ . (5) 

This result holds for both the simple periodic lattice and 
the Lieb lattice in £ e [-G, G], except right at the flat band 
energy Epp = 0 for the Lieb lattice where Eq. (1) does 
not apply. Away from this singular point, the localization 
length of the evanescent wave for both lattices is then 
given by 

f = |tm[fc]r^ = [in^cr-l- V l)] ■ (6) 

This localization length can be probed by sending a plane 
wave at energy E into the system or by introducing a point 
defect on the lattice. The latter approach was used in 
Ref. [7] where Ea in one unit cell of the Lieb lattice is 
detuned by G/10, and the resulting defect state confirms 
an evanescent tail given by Eq. (6). One can verify in the 
same way that Eq. (6) holds for the evanescent wave in 
the simple periodic lattice as well, similar in spirit to in¬ 
troducing a row defect in a 2D square lattice [31,32]. 










Importantly, we note that the localization length given 
by Eq. (6) has a minimum at the midgap [E = Eps), which 
follows the conventional understanding mentioned in 
the introduction. Although using the approach above we 
cannot evaluate the localization length at Epp on the Lieb 
lattice, it is clear that the presence of the flat band will no 
doubt make the localization length discontinuous at Epp. 
This observation can be thought as the precursor of the 
anomalous localization minimum we will discuss in the 
next section, where disorder broadens the flat band and 
smooths out the discontinuity of the localization length. 

2.2 With disorder 


that here the curvature of the localization length changes 
sign near liil = W12, where the second order derivative 
d^^[u;W)/du^ now takes a negative value. 

As we will elucidate in the next section, these anoma¬ 
lous minima are caused solely by the flat band and in¬ 
dependent of the dispersive bands on the other side of 
the band gaps. Here to highlight the anomaly of these 
localization minima, we provide another example using 
a double-sided Lieb lattice shown in Fig. 3. This lattice 
has three neighboring flat bands, and the localization 
length exhibits two minima (instead of one) inside each 
gap formed by two flat bands, which shows clearly that 
each minimum is due to the closest flat band. 


The localization length £, with disorder can be calcu¬ 
lated straightforwardly using the transfer matrix method 
(see Appendix A). Below we assume a uniform disorder W, 
with Ex j e lEx-WI2,Ex + WI2] and (Ej^j) = W^i 12. The 
results are shown in Fig. 2 as a function of the scaled en¬ 
ergy u= (E- Epp]/W and in units of the lattice constant, 
with different values of W. 

While f of the simple periodic lattice behaves simi¬ 
larly to Eq. (6) and displays a minimum atE=0, we note 
that for a given weak disorder W <G, the localization 
length f(u: W] of the Lieb lattice exhibits a minimum in 
each of the two band gaps formed by the flat band and 
its neighboring dispersive bands (i.e., ue [-G/W,0] and 
[0,G/W]). In addition, ^{u;W] barely changes with the 
disorder strength VT for | u| < 1/2, hence displaying a scal¬ 
ing behavior in terms of u. 

Importantly, the localization minima at\E\<W 12 [see 
Fig. 2(b)] are far from the respective midgap at £ = ±G/2 
when W « G, to which we refer as the anomalous lo¬ 
calization minima. To further distinguish them from the 
conventional midgap ones as shown in Fig. 2(a), we note 



Figure 3 Band structure (a) and localization length (b) of a 
ID double-sided Lieb lattice. Inset in (a): Schematic of the 
lattice structure. Solid and dashed lines indicate couplings 
G and /, respectively. Three flat band exist at Eppi = -0.5, 
Epp2 = 0 and Epp-^ = 0.5, and we have taken Ea = Ec = 0.5, 
Ep) - Ep - -0.5, Ep = 0 and G = J = 1 . The localization 
length shown in (b) is calculated with W =1/6 using the trans¬ 
fer matrix method, and the dashed lines mark the flat band 
energies. 

3 Localization length and DOS 



EA/V E/W 


To prove that the anomalous minima of the localization 
length in Fig. 2(b) are solely caused by the flat band, we 
use an integral transform of the DOS that gives the inverse 
localization length [35]. This approach utilizes the top 
right element of the Green’s matrix ^ that represents the 
probability amplitude of wave propagation from one end 
of a finite lattice with N unit cells to the other end. For the 
ID Lieb lattice, we find this matrix element to be 


Figure 2 Localization length of the lattices in Fig. 1 in the pres¬ 
ence of disorder. The vertical dashed lines show the boundaries 
of the on-site energy disorder {E e [-W/2, W/2]), and the lo¬ 
calization length is calculated with W = 0.01,0.2,0.5 and 1 
(from bottom to top). 


cofactor(£ - Hlawj ^ 

detiE-H) nfc=i(f?-Afc) 


(7) 


where H is the 3N x 3N tight-binding Hamiltonian rep¬ 
resented by Eqs. (2)-(4), {A^} are its 3N eigenvalues, and 
“cofactor" denotes tbe cofactor matrix. The localization 
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Figure 4 DOS of the Lieb lattice shown in Fig. 1 (b) with W =l 
(solid line). Its value for the flat band with W = 0.5 is shewn by 
the shadowed area, and the dotted lines mark E-±WI2. 100 
unit cells and 10,000 samples are used. 


length f (£) in this approach is given by 

lim ^ln|^i,3]vl- (8) 

N—CO N 

where the dimensionless quantity ^i, 3 n is the residual 
of ‘Si, 3 ]s[ at the eigenvalue Xk closest to E. Using Eq. (7), 
we derive a summation rule for the inverse localization 
length, i.e., 

r\E; W) = CpB (w) + ^~DB (-B: (9) 


where 


^PB^U) -f ^ 

J-2 

/ OO 

in 

-OO 


ln|u- yl Ippb[v) - l]dv, 


E — z 


pDBiz-, W)dz 


( 10 ) 

( 11 ) 


are the contributions of the flat band and dispersive 
bands, respectively. ppBM,pDBiz]W] are their respective 
DOS normalized by ppB{v)dv = I and PoBiz; W) 
dz - 1. They are related to the total DOS p{z]W) (see 
Fig. 4) hyppB{v)^3Wp[z;W){v = z/We [-1/2,1/2]) and 
pP)Biz-,W] - l.5p{z;W]{\z\ > WI2), so defined because 
there are twice as many states in the dispersive bands 
as in the flat band and we have used f^p{z; W)dz- 1. 
Equation (9) gives a good agreement with the results of 
the transfer matrix approach (see Appendix B). 

Note that we have taken (1) E - EpB + lT/2 to be the 
boundaries between the broadened flat band and dis¬ 
persive bands; and (2) ppb{v) to be independent of the 
disorder strength W, which are good approximations for 
weak disorder as exemplified by Fig. 4. The former can be 
understood since + Vk/2 is the maximum perturbation on 
each lattice site, and the latter can be understood from 
the perspective of degenerate perturbation theory: the 
flat band states are degenerate on the clean lattice, and to 


the leading order the perturbed flat band states are sim¬ 
ply their linear superpositions that are orthogonal with 
respect to the disorder potential. Clearly multiplying a 
constant to the disorder potential does not change the or¬ 
thogonality of these superpositions, and hence the energy 
of the broadened flat band states scales linearly with the 
perturbation IV, as long as the higher-order perturbations 
can be neglected. 

This perturbation argument, however, cannot be ap¬ 
plied to the localization length: even though poBiz] W) 
is zero near the flat band energy for weak disorder, 
Eq. (11) tells us that ^'j^g{E;W] here is not necessarily 
small and hence cannot be neglected. As an example, we 
show W] and its two terms in Fig. 5. We find that 

f {E; W] in fact has a larger magnitude than (u) in 
|£| < W/2, and the negligence of the former would lead 
to a negative and hence unphysical value of the localiza¬ 
tion length. Nevertheless, we note that (E; W) is al¬ 
most a constant for |£| < W/2 and weak disorder as Fig. 5 
shows. This property is due to the slowly varying kernel 
ln|£-z| inEq. (11) and the vanishing of pp,b{z]W] near 
the flat band energy for weak disorder (see Fig. 4). Hence 
W], which represents the coupling to the disper¬ 
sive bands, does not play a role in the non-monotonic en¬ 
ergy dependence of f (£; W] near EpB, and the anomalous 
minimum of f (E; W] on each side of EpB is solely due to 
the flat band itself: as f (E) in Fig. 5 clearly shows, these 
minima correspond to the maxima of (E), and the lat¬ 
ter can be captured by a crude (and linear) approximation 
for the DOS of the flat band, i.e., ppb{u) 2-4|u|(|w| < 
1/2). Once substituted into Eq. (10), it gives the following 
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Figure 5 Inverse localization length of the Lieb lattice shown 
in Fig. 1(b) with W = 0.5 (thick solid line). The dash-dotted 
and dashed lines show the flat band and dispersive band term 
in Eq. (9), respectively. The dotted line shows the analytical 
approximation (12) in juj < 1/2 (delimited by the vertical lines), 
and the arrows mark its peak positions at juj = 0.46, which 
agree well with those of the numerical results. 


















analytical expression 




-h |u|ln 

2 


i + |m| 
i-|M| 


+ 2u^ In 



( 12 ) 


for I m| < 1 /2, which agrees reasonably well with its value 
obtained numerically (see Fig. 5), including its peak posi¬ 
tions (and hence also the positions of the minima of the 
localization length) near \ u\ = 0.46. 

We further note that W) is also insensitive to a 

weak disorder W for a fixed E (in addition to E for a given 
W mentioned above) near the flat band, again due to the 
slowly varying kernel InlU- z| in Eq. (11). It smears out 
the change of posiE; W] when W changes to W', i.e., 


^-'g(E;W]-^^]^[E;W']^ f In 

Jdb 


E-z 


dz 

T 


• f [pDB (z; W) - pdb iz; W']]dz^0, (13) 

Jdb 


when the energy E is far from the dispersive bands. The 
subscript “db " of the integrations denotes integrating over 
the bandwidth F of the dispersive bands, and the above 
difference vanishes due to the second integral and the 
same normalization of pp)b(E< V/),poBiE; W']. The van¬ 
ishing difference above explains why the different local¬ 
ization curves in Fig. 2(b) overlap in | u| < 1/2, which leads 
to a scaling behavior 


aE-,W)^au) (|u|<l/2) 


(14) 


for weak disorder. We note that the W-independence of 
the localization length right at £ = Eps (and u = 0) found 
previously [18,22,23] is a special case of Eq. (14). 


4 Flat band states at the localization 
minimum 

In the previous section we have established that the 
anomalous minimum of the localization length on each 
side of EpB is solely due to the DOS of the flat band and 
independent of the opposing band. In this section we dis¬ 
cuss the spatial profile of the states at these localization 
minima, which we found to be the compact states that 
exist even in the absence of disorder [2,18]. 

We have noted that E = EpB ±WI2 delimit the non¬ 
zero range of Ppb{u] when discussing Fig. 4. This obser¬ 
vation is in fact quite counterintuitive but also instructive: 
the maximum energy shift of each lattice site is W/2, and 
in order to shift the energy of a state spanning several lat¬ 
tice sites by this amount, all these lattice sites then need to 
have the same disorder potential W/2. The probability of 
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Figure 6 Top: Example of on-site disorder in a Lieb lattice with 
150 sites (50 unit cells) and W = 0.5. The value of the disorder 
on each lattice site is illustrated using the color scale on the 
right. Bottom: Spatial profiles of three representative states in 
this example, including a compact flat band state near u = 0.41 
(dark solid line), a more extended flat band state near u = 0.08 
(red dashed line), and a dispersive band state (grey solid line). 
Inset: Schematic of the compact state overlaid with the lattice 
structure. 


such a rare event is negligible for modes in the dispersive 
band, which have a long localization length (see the grey 
line in Fig. 6). In the flat band however, there are compact 
states that span only two unit cells in the absence of disor¬ 
der [2,18], and they are “dark" on the B sites where their 
wave functions vanish due to the cancellation of hopping 
from neighboring A and C sites (see Fig. 6; inset). Hence 
these compact states are insensitive to the disorder on 
the B sites and only require three similar on-site disorder 
potentials on two consecutive A sites and the C site in 
between to have a nearly unaltered compact state. One 
example of a compact state at u = 0.41 is shown in Fig. 6, 
and we note that for a uniform disorder, the probability 
of having such compact states near | u| = 1/2 is the same 
as that near \ u\ =0, which leads to a long tail of pppiu) 
(up to |u| = 1/2) when compared with that of poBiE; W] 
inside the band gaps (see Fig. 4). 

These compact states are weakly coupled to the rest 
of the disordered lattice and lead to the minima of ^(u) 
near | u| = 1 /2: the flat band modes other than these com¬ 
pact states are more extended along the lattice, and they 
sample more on-site disorders which tend to cancel each 
other (see the flat band mode at u = 0.08 shown by the 
dashed line in Fig. 6). Therefore, the energies of these 
more extended flat band modes are closer to the undis¬ 
turbed flat band energy, which leads to the increase of the 
localization length as \ u\ approaches 0. In |u| > 1/2, the 
flat band states do not exist and ^pg(u; W) approaches 




















Figure 7 Band structure (a) and localization length (b) of a 
1D Lieb lattice with a bent band in the middle. In (a) all the 
parameters are the same as in Fig. 1(b) except for Ea = 0.1. 
The center band has maxfil) = 0.08 and min(£) = 0, and the 
dashed line shows the flat band when Ea = 0. Its localization 
length in (b) is calculated with W = 0.2 (solid line) and 0.8 
(dash-dotted line), (c) and (d) Typical modes on the high and 
low energy sides of the bent band with W = 0.2, respectively, 
(e) and (f) Same as (c) and (d) but with W = 0.8. 


zero (see Fig. 5), and as a result, the localization length 
gradually increases to its value given hj ^db(.E;W] due to 
the dispersive hands. 

To further emphasize the role of the compact states in 
forming the anomalous minima of the localization length 
near we note that these states persist even when the 
flat band is slightly bent, for example, by introducing a 
small detuning e > 0 to the A sites [see Fig. 7(a)]. Note 
however, they now only appear on the high energy side 
of the bent band when W12 <e [see Figs. 7(c) and (d)], 
and hence the anomalous minimum of the localization 
length should only appear on the high energy side also, 
as we verify in Fig. 7(b): the minimum of the localization 
length on the low energy side of the bent band, in con¬ 
trast, follows the conventional scenario due to the lack of 
the compact states, and the curvature of the localization 
length does not change its sign inside this band gap. As 
the disorder strength increases, i.e., when W/2 becomes 
large when compared with e, the compact states appear 
on both sides of the bent band [see Figs. 7(e) and (f)], and 
the anomalous minimum of ^(£; W] now reappears on 
the low energy side of the bent band [see Fig. 7(b)]. 


5 Conclusion and Discussion 

In summary, we have shown that anomalous localization 
minima near an isolated flat band is a generic feature 
caused by its sharp DOS. In the previous instance of these 
minima found in a photonic crystal waveguide, they exist 
for evanescent waves in a defect-free system, where the 
flat band is due to a Van Flove saddle point singularity. 
Here we have shown that they exist in the presence of 
disorder in a Lieb lattice, where the flat band is the result 
of canceled quantum tunneling from A, C sites to the 
neighboring B site. 

As a result of this generic property, when a more exotic 
band gap is formed by two flat bands, we find two such lo¬ 
calization minima inside this band gap, which highlights 
their anomaly in view of the conventional midgap min¬ 
imum of the localization length exemplified in Fig. 2(a). 
These anomalous minima are direct manifestations of 
the compact states in the flat band, which persist even 
when the flat band is slightly bent. In the examples we 
have shown the on-site energy is subjected to a uniform 
disorder, and we have checked that these features hold for 
other types of disorder as well, such as a Gaussian disorder 
(see Appendix C). 

As mentioned previously, the tight-binding model 
given by Eqs. (2)-(4) can be applied to different types 
of particles as long as the noninteracting picture ap¬ 
plies. Therefore, although the variation of the localization 
length near Efb is on the order of the lattice constant and 
hence difficult to observe for electrons, it can be probed in 
synthesized systems as numerically shown in Refs. [31,32], 
using a quasi-ID dielectric photonic crystal waveguide. 
The candidates to observe these anomalous localization 
minima on a Lieb lattice include cold atoms in optical 
lattices and exciton-polariton condensates in miropillar 
cavities [17,18], as long as the interaction between the 
underlying particles is kept small. These systems have a 
lattice constant on the order of microns [17,18] or tens 
of microns [4-6], and hence they are ideal for verifying 
the anomalous minima of the localization length (see the 
experimental technique used in Ref. 18, for example). 

We have also shown that the localization length dis¬ 
plays a scaling behavior close to Efb, which holds even 
when the flat band is shifted toward one dispersive band 
(see Appendix B). Using the top right element of the 
Green’s matrix, we have derived an analytical relation be¬ 
tween the DOS and the localization length for the ID Lieb 
lattice. It takes the form of a summation for the inverse 
localization length, with one flat band term and one dis¬ 
persive band term, and we have used it to explain the 
localization properties mentioned above. 
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Finally, we note that while localization in ID disor¬ 
dered systems is qualitatively different from 2D and three- 
dimensional systems in general [36], the anomalous min¬ 
ima of the localization length near an isolated flat band 
are tied to the compact states which exist in higher dimen¬ 
sions as well [2,23]. Hence it is likely that these anomalous 
minima may exist beyond ID, even though they have not 
been found in 2D systems with a flat band touching a 
dispersive band [20]. 
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A Transfer matrix calculation 


The transfer matrix method is a convenient way to cal¬ 
culate the localization length in ID [37]. It relies on the 
following relationship between the values of the wave 
function on three consecutive lattice sites of the same 
type: 


Xj+i - djXj + ej-iXj^i. 


(15) 


The transfer matrix tj for neighboring unit cells is then 
defined by 


Y 

_ 

dj Bj-i 

(xj] 

= ■*■/ 


\xjj 


il 0 I 

U-iJ 

J 

k-J 


and the total transfer matrix is given by 
JV-l 

T=Y\^j’ (17) 

i=2 


satisfying 

' Xn 
Xn-1) 



(18) 


The exponent of its largest eigenvalue Amax. after ensem¬ 
ble average (denoted by {■» and divided by N, gives the 
inverse localization length in the asymptotic limit 
N —>■ oo: 

lim 4(ln(Amax)>- (19) 

JV— CO N 


The recursive relation (15) can be easily found for A, 
B,orC sites of the ID Lieb lattice shown in Fig. 1 (b) of the 
main text. Take B sites for example, by expressing Aj, Cj 
in terms of Bj and Bj+i and plugging them into Eq. (3) in 
the main text, we find 


fj ~ - - Xj ~ Xj-1’ Xj - T > 

^Aj ^Cj 


( 20 ) 

( 21 ) 


Figure 8 Localization length in a 1D Lieb lattice atW= 0.05. 
The solid line and open circles are calculated using the recur¬ 
sive relation (20) and (22), respectively. The other parameters 
are the same as in Fig. 1 (b) in the main text. 


or equivalently dj = fjIxj, Sj-i = -Xj-i/Xj in form of 
Eq. (15). This is the option we have chosen in the main 
text. 

Similarly, we find the following recursive relation for C 
sites 


fjCj - Xj+iCj+i + XjCj-i, (22) 

f 

f] = Xcj-Xj-Xj+i’ Xj = -(23) 

which gives the same localization length as the first choice 
above (see Eig. 8). 


B Compare localization length obtained 
using different methods 

The results obtained from the transfer matrix and the ana¬ 
lytical expression (9) in the main text agree well, which we 
illustrate in Eig. 9. In Fig. 9(a) the flat band sits at Epg = 0 
as in Fig. 1 (a) of the main text, and in Fig. 9(b) we shift the 
flat band to£i7B = 0.1 = e without changing the dispersive 
bands. The latter is achieved by choosing Ea = Ec = e, 
Eb - -e and G —>^ \/G^ We note that the scaling be¬ 
havior given by Eq. (14) in the main text holds near Epg 
even though Epg is shifted. This result is consistent with 
the proof of Eq. (14) we have given, and all it requires is 
that the dispersive band contribution W) is insen¬ 

sitive to W, which is satisfied when the shifted Epp is still 
far from the dispersive band edges with weak disorder. We 
do note that f (E; W) now shows a slight asymmetry about 
E - EpB ^0{u = 0), which comes from the contribution of 



















the two dispersive bands, which are no longer symmetric 
with respect to the shifted flat band. 
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Figure 9 (a) Localization length of the Lieb lattice shown in 
Fig. 2(b) of the main text (solid lines), calculated with W = 
0.5 (bottom curve) and 1 (top curve) using the transfer matrix 
method. The dashed lines show the results using Eq. (9) in 
the main text and numerically calculated DOS (see Fig. 4 for 
example), (b) Same as (a) but with the flat band shifted to 
= 0 . 1 . 


C Localization length with a Gaussian 
disorder 

In the conclusion section of the main text, we mentioned 
that the anomalous minima of the localization length near 
an isolated flat band occur for other types of disorder 
as well. Here we show two examples using a Gaussian 
disorder, one for the Lieb lattice in Fig. 1 (b) and one for the 



Figure 10 Anomalous minima of the localization length with a 
Gaussian disorder, for (a) the Lieb lattiee in Fig. 1(b) and (b) 
the double-sided Lieb lattice in Fig. 3. We express the variance 
of the on-site disorder in the same form as that for a uniform 
disorder, i.e., (E^,) = W^/12, where W = 0.1 in (a) and 1/6 in 
(b). 


double-sided Lieb lattice in Fig. 3. The results are shown 
in Fig. 10, where two and six such minima are clearly seen 
in these two cases, respectively. 

Acknowledgements. The author thank Hakan Tiireci, 
Vadim Oganesyan, Matteo Biondi, Sebastian Schmidt and 
Bo Zhen for stimulating discussions. This project is sup¬ 
ported by NSF under Grant No. DMR-1506987. 

Key words, flatband, localization length, Lieb lattice. 


References 

[1] V. Apaja, M. Hyrkas, and M. Manninen, Phys. Rev. A 
82, 041402(R) (2010). 

[2] M. Hyrkas, V. Apaja, and M. Manninen, Phys. Rev. A 
87, 023614 (2013). 

[3] M. C. Rechtsman, ]. M. Zeuner, A. Ttinnermann, S. 
Nolle, M. Segev, and A. Szameit, Nature Photon. 7, 
153 (2013). 

[4] R. A. Vicencio, C. CantiUano, L. Morales-Inostroza, B. 
Real, C. Mejia-Cortes, S. Weimann, A. Szameit, and M. 
1. Molina, Phys. Rev. Lett. 114, 245503 (2015). 

[5] S. Mukherjee, A. Spracklen, D. Choudhury, N. Gold¬ 
man, P. Ohberg, E. Andersson, and R. R. Thomson, 
Phys. Rev. Lett. 114, 245504 (2015). 

[6] D. Guzman-Silva, C. Mejia-Cortes, M. A. Bandres, 

M. C. Rechtsman, S. Weimann, S. Nolle, M. Segev, A. 
Szameit, and R. A. Vicencio New). Phys. 16, 063061 
(2014). 

[7] L. Ge, Phys. Rev. A 92, 052103 (2015). 

[8] C. L. Kane and E. ]. Mele, Phys. Rev. Lett. 78, 1932 
(1997). 

[9] F. Guinea, M. 1. Katsnelson, and A. K. Geim, Nature 
Phys. 6, 30 (2010). 

[10] A. Simon, Angew. Chem. 109, 1873 (1997). 

[11] S. Deng, A. Simon, and). Kohler, Angew. Chem. 110, 
664 (1998). 

[12] S. Deng, A. Simon, and 1. Kohler, L Solid State Chem. 
176,412 (2003). 

[13] M. Imada and M. Kohno, Phys. Rev. Lett. 84, 
143(2000). 

[14] E. Tang, 1-W. Mei, and X-G. Wen, Phys. Rev. Lett. 106, 
236802 (2011). 

[15] T. Neupert, L. Santos, C. Chamon, and C. Mudry, 
Phys. Rev. Lett. 106, 236804 (2011). 

[16] S. Yang, Z.-C. Gu, K. Sun, and S. Das Sarma, Phys. Rev. 
B86, 241112(R) (2012). 

[17] T. Jacqmin, 1. Carusotto, 1. Sagnes, M. Abbarchi, D. D. 
Solnyshkov, G. Malpuech, E. Galopin, A. Lemaitre, 

1. Bloch, and A. Amo, Phys. Rev. Lett. 112, 116402 
(2014). 

[18] F. Baboux, L. Ge, T. Jacqmin, M. Biondi, E. Galopin, 

A. Lemaitre, L. Le Gratiet, 1. Sagnes, S. Schmidt, H. 

E. Tiireci, A. Amo, and J. Bloch, Phys. Rev. Lett. 116, 
066402 (2016). 





















[19] M. Goda, S. Nishino, and H. Matsuda, Phys. Rev. Lett. 
96, 126401 (2006). 

[20] J. T. Chalker, T. S. Pickles, and P. Shukla, Phys. Rev. B 
82, 104209 (2010). 

[21] J. D. Bodyfelt, D. Leykam, C. Danieli, X. Yu, and S. 
Flach, Phys. Rev. Lett. 113, 236403 (2014). 

[22] D. Leykam, S. Flach, O. Bahat-Treidel, and A. S. Desy- 
atnikov, Phys. Rev. B 88, 224203 (2013). 

[23] S. Flach, D. Leykam, J. D. Bodyfelt, P. Matthies, and A. 
S. Desyatnikov, Europhys. Lett. 105, 30001 (2014). 

[24] J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. 
D. Meade, Photonic Crystals: Modeling the Flow of 
Light, 2 ed. (Princeton University Press, New Jersey, 
2008). 

[25] P W. Anderson, Phys. Rev. 109, 1492 (1958). 

[26] F. Urhach, Phys. Rev. 92, 1324 (1953). 

[27] G. D. Cody, in: Hydrogenated Amorphous Silicon 
Optical Properties, edited hy B. J. Pankove, Semicon¬ 
ductors and Semimetals, Vol. 21B (Academic, New 
York, 1984). 

[28] N. F. Mott and E. A. Davis, Electronic Processes in 
Non-Crystalline Materials, 2 ed. (Oxford University 
Press, New York, 2012). 

[29] M. H. Cohen, E. N. Economou, and C. M. Soukoulis, 
Phys. Rev. B 32, 8268 (1985). 

[30] E. N. Economou, C. M. Soukoulis, M. H. Cohen, and 
A. D. Zdetsis, Phys. Rev. B 31, 6172 (1985). 

[31] Z.-Y. Li and K.-M. Ho, Phys. Rev. Lett. 92, 063904 
(2004). 

[32] M. Ibanescu, E. J. Reed, and J. D. Joannopoulos, Phys. 
Rev. Lett. 96, 033904 (2006). 

[33] A. Lagendijk, B. van Tiggelen, and D. S. Wiersma, 
Phys. Today 62, 24 (2009). 

[34] T. Schwartz, G. Bartall, S. Fishman, and M. Segev, 
Nature 446, 52 (2007). 

[35] D. C. Herbert and R. Jones, J. Phys. C: Solid St. Phys. 4, 
1145 (1971). 

[36] P. A. Lee and T. V Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

[37] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 
1469-1564 (1993). 


